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We study quantum dynamics of the rotationally driven Dicke model where the collective spin is 
rotated around the z axis with a finite velocity. In the absence of the rotating wave approximation we 
observe that for several physically relevant initial states the position of the quantum critical point is 
shifted by the amount given by the applied rotation velocity. This allows us to probe the quantum 
criticality "from a distance" in parameter space without actual crossing of the quantum critical 
surface but instead by encircling it in the parameter space. This may provide a useful experimental 
hint since the quantum state is not destroyed by this protocol. Moreover, for the coherent initial 
state we observe an interesting non-equilibrium reentrant phenomenon of quantum critical behavior 
as a function of the driving velocity and construct a non-equilibrium phase diagram of the driven 
model. 



PACS numbers: 

I. INTRODUCTION 

The Dicke model pQ, [2] was introduced to model the 
interaction of a collection of N two-level atoms with a 
single-mode radiation (bosonic) field through the dipole 
coupling. Under the condition that the ensemble of two- 
level atoms is confined in a region of space which is much 
smaller than the wavelength of the single-mode radiation 
field, the system of the two-level emitters behaves as a 
collective large spin. The model became a paradigm for 
collective behavior of matter interacting with light and is 
believed to have a phase transition at finite and zero tem- 
peratures [3], [4], [5]. A super-radiant transition is char- 
acterized by the expectation value of the photon num- 
ber operator which distinguishes between normal phase 
(where it is zero) and a super-radiant phase (where it 
is nonzero). The issue of this phase transition has been 
debated [B] in length in the literature by various authors 
using different methods. The conclusion is that although 
is it difficult to realize the phase transition in a real two- 
level system for which the model has been introduced, it 
is possible to engineer it in some cavity QED setups. For 
a recent review of these issues see, e.g., [7], [H]- 

With a recent progress in quantum optics and cold 
atoms it became possible to realize the Dicke model in a 
cavity QED systems and observe a super-radiant phase 
transition and other collective behaviors [9^. This ad- 
vance became possible by achieving considerable degree 
of tunebility of parameters of the system with high pre- 
cision. Moreover parameters can be controlled in real 
time. This stimulates studies of non-equilibrium effects 
in this class of systems. An interplay between collec- 
tive and non-equilibrium effects is pronounced close to 
the quantum phase transition. Several recent theoreti- 
cal works on driven non-equilibrium Dicke model, |10j . 
[TT] , [12] and [13] have focused on various effects of time- 
depending driving of the interaction strength and/or de- 
tuning. Many of these non-equilibrium effects are related 
to the presence of a quantum critical point in the equi- 
librium model. 



Here we propose to probe the quantum phase transi- 
tion appearing in the Dicke model in an indirect way. 
We employ a certain degree of tunebility of the system 
to drive (rotate) the collective spin by a certain time- 
dependent angle. This makes it possible to probe the 
quantum criticality "from a distance" in parameter space 
without crossing the equilibrium quantum critical point. 
Apart from this we observe new, reentrant quantum crit- 
ical behavior in such non-equilibrium protocol for certain 
initial state. This behavior can be traced back to a mech- 
anism of competition between geometric and dynamical 
phases in the non-equilibrium dynamics |14| . 

In addition to the previous studies, here we look into 
the quantum dynamics of the rotationally driven Dicke 
model when the change of the rotation angle is linear 
in time. We distinguish between two types of protocols 
(see Fig. [T]): in one case we encircle the critical surface 
in the parameter space from the outer region (contour 
Ci) while in the second case we drive the system in such 
a way that the circle C2 stays inside the region bounded 
by the critical surface. These regions are distinguished 
by the presence (region Ci) or absence (region C2) of the 
geometric phase contributions to the total phase factor of 
each instantaneous eigenstate. To characterize the non- 
equilibrium phase diagram we look into two types of ob- 
servables: (i) the mean photon number and (ii) a topo- 
logical number associated with the expectation value of 
the parity operator. 

Moreover we consider several different initial states: 
(i) atom-field coherent states (coherent state for photons 
and pseudo-spin subsystem), (ii) Fock state and (iii) the 
ground state of the system. Our quasi-classical approach 
which is based on the coherent state representation is 
supplemented by a numerical procedure. For initial co- 
herent and ground states a non-equilibrium phase dia- 
gram has an additional reentrant phase inside the super- 
radiant phase. For the aims of completeness and com- 
parison with a driven system, we also look into dynam- 
ics of the undriven Dicke model starting from those ini- 
tial states. We study the system both for a finite num- 
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Figure 1: Equilibrium quantum phase diagram of the rota- 
tionally driven Dicke model as a function of atomic frequency 
ujo and atom-photon coupling strength A. Paths Ci and C2 are 
fixed by the time dependence of the driving parameter (j){t) 
(rotation angle around z-direction in spin space). Here the 
bosonic frequency is set to a; = 1. 



ber of two-level systems and in the thermodynamic limit 
(TDL). We found in particular that dynamics starting 
from the coherent state coincides with the one for the 
initial ground state, both in the driven and undriven 
cases. Dynamics from the pure Fock state is featureless 
in the thermodynamic limit: one need some fluctuations 
in the initial state to make it evolve, at least in the quasi- 
classical approach. This fluctuations are introduced by 
allowing small fluctuations in the bosonic and the spin 
sectors; we call it "nearly-Fock" state. Our results for 
the non-equilibrium phase diagram of the driven Dicke 
model are summarized in Fig. [8] 

We note that here we study the Dicke model beyond 
the rotating wave approximation which means that it is 
not integrable. Dynamics of integrable variant of the 
Dicke model starting from certain initial states (when the 
rotating wave approximation is made) has been studied 
before in [13] and more recently in [TB] , [17] . We however 
found some features (like e.g. formation of soliton-like 
trains in the thermodynamic limit) which also hold in 
the non-integrable case. 

First, we overview some basic facts about Dicke model 
(Section II), while in the Section III we establish quasi- 
classical equations for the dynamically driven model. 
Section IV presents various results for several initial con- 
ditions and comparison between dynamics of the driven 
and undriven models. 



II. THE DICKE MODEL 

Here we consider the Dicke model without the rotating 
wave approximation. The quantum Hamiltonian of the 



Dicke model reads 
Hd = t^oJ^ + i^a^a - 



A 

N' 



(a) + a) (j 



(1) 



where ujq is the level-splitting of a two-level atom and 
the collection of those N identical two-level atoms is de- 



scribed by the collective spin operators — J2iLi l'^- 
Here , /.t = z, ± are the usual Pauli matrices which de- 
scribe a two-level emitter. The collective spin operators 
satisfy the su(2) Lie algebra: 



[j+,j_] = 2j„ [j„j±]=±J±. 



(2) 



The Hilbert space of the su(2) Lie algebra is spanned by 
the Dicke states |j, m), with m = —j, — j -I- 1, . . . , j — 1, j 
and j being integer or half integer. The Dicke states 
I j, m) are the eigenstates of and Jz with eigenvalues 
+ s-iid TO, respectively. As in |18j we fix j to its max- 
imal value j = N/2. Therefore the collection of two-level 
atoms can be interpreted as a large "spin" (pseudo-spin) 
vector of magnitude j = N/2. The single-mode radia- 
tion field with frequency cj is described by the bosonic 
creation and annihilation operators a, a^ satisfying the 
usual commutation relation [a, a^] = 1. The set of opera- 
tors {a^a, a^, a, 1} forms the Heisenberg-Weyl Lie algebra 
h4. The Hilbert space of this algebra is spanned by the 
Fock states |n) with n = 0, 1, 2, 3, . . ., defined by 



a^lTi) 



1 



1) 



an) = \/n\n — 1). 



(3) 



Finally, the term -;=,(§,' + a)(J+ + J_) describes the 



atom-field interaction within the dipole approximation 
where the coupling strength is A and the factor stems 
from the fact that the original dipole coupling strength 
is proportional to 

The symmetries of the Dicke model are found by ap- 
plying a unitary transformation of the form 



U((y5o) = exp icpolSl 



(4) 



to the Dicke Hamiltonian Hp. Here N = a)a. + + j 
is the excitation number operator, the operator which 
counts the number of excited quanta in the system and 
ipo is some rotation angle. Applying this transformation 
we observe that 



UaU^ 



(5) 



and therefore, the transformed Dicke Hamiltonian reads 

UHdU^ - uJoh 
A 



cja^a 



+ 



/2J 



(a^J+e^*'^" + a^J_ + aj+ + aJ^e-^*'^") . (6) 



The Dicke model remains invariant under this transfor- 
mation only if ifQ = TT. Hence, the associated symmetry 
of the Dicke model without the rotating wave approxi- 
mation is the parity, since 



V{ipo — tt) = exp iirN 



H. 



(7) 
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This has the interpretation as a Z2 symmetry because it 
has the eigenvalues e"("+'"+J) = ±1. Thus the Hilbert 
space breaks up into two non-interacting subspaces, one 
with an even excitation number n -\- m + j and one with 
an odd one. Therefore, the Dicke model without the 
rotating wave approximation possesses a Z2 symmetry. 

The Dicke model ([ij has a quantum phase transition 
in the thermodynamic limit N ^ 00 {j ^ 00) at the 

quantum critical point Ac — . Quantum critical 

points (criticalities) are points in the underlying param- 
eter space (cjo,cj, A) at which the excitation energy van- 
ishes. Hence, the ground state of the system changes 
qualitatively from a normal to a super-radiant one. In 
the super-radiant phase the ground state is a macroscop- 
ically excited collective state that possesses the poten- 
tial to super-radiate. A phase transition is associated to 
spontaneous breakdown of the Z2 symmetry described 
above. The super-radiant emission is a collective effect 
involving all the N two-level atoms in the sample. In this 
case the decay rate is F ^ ^^70 instead of A^7o, the rate 
for independent atom emission. Here 70 is the decay rate 
of an excited two-level atom. 

In this paper we study a rotationally driven version of 
the Dicke model [T5j, [5^, which results from applying a 
rotation of angle around the z axis to each two-level 
atom. Explicitly it has the form of time-dependent uni- 
tary rotation applied to z axis, 

R,it)^exp -i(j){t)h ■ (8) 

The resulting time dependent Hamiltonian HpiD(t) — 
Rt(t)HD(wo,a;,A)R^(t) reads 

+ (at + a) (e'^(*) j+ + e^*^(*)j_) . (9) 

Hence, the applied time dependent rotation (t) can be 
interpreted as a rotation of the pseudo-spin vector J of 
length j = N/2 around the z axis by an angle (j>(t). 

Since the rotated Dicke Hamiltonian HrdIO is ob- 
tained by applying a unitary transformation, namely by 
the rotation R^ (t) , the eigenvalues (eigenenergies) of the 
rotated Hamiltonian remain the same. Therefore, the 
quantum criticalities of the rotated Dicke Hamiltonian 
can be obtained by applying a rotation around the z axis 
to the time independent Dicke Hamiltonian. We use this 
fact and the results of Ref. 18 where the time indepen- 
dent Dicke Hamiltonian is diagonalized in the thermody- 
namic limit. In the "normal" phase (NP) the excitation 
energy reads 



(NP) _ 




< 



in which the reality condition A ^ ^ 
fulfilled. Similarly, in the "super-radiant' 



(10) 

needs to be 
phase (SRP) 



the excitation energy is 
(SRP) _ 1 / ISA'^ 



(11) 



which is real under the condition that A > , where 

f{uj, A) = (^Ifr- — -1-4^2^2^ Hence, the critical cou- 
pling at which the quantum phase transition occurs is 
given by Ac = . Using this relation we plot the 

equilibrium phase diagram for uj — 1. We parametrize it 
by A cos (j) for the x axis, and A sin (p for the y axis and 
ujQ on the z axis. This is shown in Fig. [l] for the purpose 
of illustration of our setup. 



III. COHERENT STATE REPRESENTATION - 
CLASSICAL ANALOGUE - TIME DEPENDENT 
MEAN FIELD THEORY 

For concreteness of our protocol we will study quantum 
dynamics governed by the Hamiltonian ^ when (j){t) is 
linear in time 



(f){t) = Sc/, t, 



(12) 



where 6^ is the driving velocity and the time interval 
is < t < tf. This driving protocol generates circular 
paths in parameter space (see Fig. [T]). For tf — 27r/(50 
the path corresponds to one closed circle. First we will 
be interested in the mean photon number (a^a) for the 
closed circular paths Ci and C2 



stsi 



(13) 



and for different initial quantum states liplO)). Time evo- 
lution is governed by the time dependent Schrodinger 
equation idt\ip(t)) — HRD(i)|'0(O)- Unfortunately, the 
Hamiltonian ([9| HRD(i) is not a linear function of the 
generators of the dynamical group, here SU(2) and 
H4. Therefore, an exact solution of the time dependent 
Schrodinger equation is difficult to obtain and we have 
to rely on some approximate methods which may become 
exact in certain limits (e.g. in the thermodynamic limit 
of j 00). For this reason we will be interested in the 
dynamics in the large spin limit (j — ^ 00) and also since 
the actual quantum phase transition only occurs in the 
thermodynamic limit. In this respect the use of the basis 
of coherent states is optimal because of the semi-classical 
nature of this basis. Further, the use of the coherent 
states allows us to write some time dependent mean field 
approximation for the dynamics which in the thermo- 
dynamic limit (semi-classical limit) becomes exact. In 
particular, the spin coherent states [H] (SU(2) coherent 
states) are defined as 



m) = 



1 



exp 



(cWJ+)b-,~j) (14) 
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where is the lowest eigenstate of Jz, Jz|j, — i) = 

The parameter C(t) G C results from the 
parametrization of the two-dimensional sphere by the 
stereographic projection 



C = tan-e*'^. 



(15) 



Matrix elements of the spin operators in this basis are 
given by 



(C|J-|C)-((C|J+IC))*, 
(C|J.|C) = C(C|J+IC>-j = -j 



1-CC 



1 + CC* 

Similarly, for the bosonic part of the system we have 



(16) 
(17) 
(18) 



\a{t)) = cxp (^-^a{t)a* {t)^ exp {a{t)k^) |0), (19) 

where |0) is the field vacuum state a|0) = and a{t) £ C 
a complex number. It is also convenient to introduce the 
following parametrization 



qi{t)+ipi{t) 



a{t) = 4?('Z2(i)+*P2(t)), 



(20) 
(21) 



where {qi,pi,q2,p2) describe the phase space of the sys- 
tem under consideration and the indices 1 and 2 label the 
spin and the field subsystem, respectively. 

Following the standard procedure of [3T], [22] we can 
obtain time dependent mean field equations of motions. 
Namely, defining the classical Hamiltonian as 

(a|(C|HRD(t)|C>|a> =^ci(0 

= ^^0 (<Z? +pI- 2j) + ^cj (ql + pI) + 



_j_2^^4^J|i+pi) (g^cos0(t)+pisin,/)(i))g2. 

(22) 

we consider a minimization problem for the action inte- 
gral given by 



S = 



(23) 



where the corresponding Lagrangian is given by 

£((?!, gi, 92,92,*) = (a|(C|«9t-HRD|C)|a)- (24) 

The minimization leads to the Euler-Lagrange equations 
of motions 



d rdc 



dt \ dqi J dq. 



dC 



= 0, 



(25) 



tion 



from which we obtain the Hamiltonian equations of mo- 
— —pi and ^^^^ = qi taking the following form 

{cos(f>{t)qi +sm(f){t)pi)piq2 



dq 

qi = uopi ~ 2A 



2A >'-(g.+^^) sin0(O,2, 



(26) 



Pi = -t^o^i + 2A 



qi {cos(j){t)qi +sm(l){t)pi)q2 



4j 



92 = UJP2, 
P2 



(27) 
(28) 



ujq2-2xj^ (^LtPfl (cos</)(i)9i+sin0(i)pi) 



4j 



(29) 



First, let us determine the fixed points of these 
equations of motions in the classical phase space 
{qi,Pi,q2,P2)- Since, these equations are explicitly time 
dependent we need to apply a canonical transformation 
(91, Pi, 92, P2) -> (Qi, A, 92,^2) that removes the explicit 
time dependence and then find out the fixed points. The 
first and the simplest fixed point is 



91, ci = 0,pi,ci = 0,92,ci = 0,P2,ci = 0. 

This fixed point is only stable if 



A < 



In addition, there are two fixed "circles" given by 



and 



9l,C2 
Pl,Cl 

P2,C2 
Pl,C3 = 

I2,cs = 
P2,cs = 



4A2 



I2j 



1 j sm<^(i). 



2A 

CO 

0, 



i 1- 



4A2 



4A2 



I2j 
2A 

UJ 



1 j sm^it), 



\ 



4A2 



(30) 
(31) 

(32) 

(33) 

(34) 
(35) 

(36) 
(37) 

(38) 
(39) 
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which are only real and stable if 



A > 



(40) 



We plot the phase portrait of the system in Fig. [2j 



Wo=l, cj=i, A = l, j = 1000, (5^=1 




Figure 2: Plot of the phase portrait {qi{t),pi{t)) for different 
initial conditions. 



displacement parameters a(t) and f3{t): 



ice — Loa — Ai 



2j 



ict - uja* - \J ^-^ ^^^* (/3*e*^ + ^e"*^) = 0, 



(46) 



(47) 



i/3 - LUof3 - A^/ ^-^ J'^* {a + a*) e'*+ 



2j 



1 



2 V2j - 13(3* 



(a + a*) {P*e'^ + ^e-"'') 13 ^ 0, 

(48) 



i/3* -ujo/3* - A^ 
1 



l2j 



2j 



(a + a*)e-^'^+ 



+ -A-L^^ " (a + a*)(p*e'''' + I3e~"*') (3* = 0, 

(49) 

which are equivalent to the system of equations obtained 
by the mean field approximation (26)-(29) if we use the 
following identification 

{qi + ipi) , a = -4^ ('72 + ip2) ■ (50) 



It is instructive to note that the classical equations 
above can also be derived following a different proce- 
dure. Using the approach of Ref. [TB] one can introduce a 
Holstein-Primakoff representation for the spin variables 
as 



2j 



-tjtb. 



j- = V2jyi-^btb b, 

1 = 9h - 1. 



(41) 

(42) 
(43) 



Next we apply a time-dependent displacement to both 
the bosonic field operators a and the Holstein-Primakoff 
bosons b, 

f)^{a{t)) = exp {a{t)&^ - a*{t)a) , (44) 
i>^il3{t)) = exp (/3(t)bt - (3*{t)h) . (45) 



This transformation takes into account the fact that the 
bosonic field operators and the atomic ensemble may ac- 
quire a macroscopic occupations. In the thermodynamic 
limit (j — >■ oo) we find the following equations for the 



IV. RESULTS FOR THE MEAN PHOTON 
NUMBER AND THE PARITY 



In this Section we study the time evolution of the 
scaled mean photon number and the parity for the ro- 
tated Dicke model. We explore several initial states and 
compare the results between driven and undriven situ- 
ations. Namely, we consider several physically-relevant 
initial states: (i) coherent states |a)|C), (ii) Fock state 
\n = 0)|j, m — —j) and (iii) the ground state |GS). Be- 
fore going into detailed studies we introduce important 
observables and review the details of numerical proce- 
dure. 



A. General setup 

1. Mean photon number 

The scaled mean photon number for the rotationally 
driven Dicke model is defined by 



1 



5ts\ - 



1 



a'a) = -{ij{tf)\sJa\iP{tf)) 



= -mo)\U^{tf)a^aU{tjmto)}, (51) 
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where the evolution operator is the time ordered expo- 
nent 



Kit) = Tcxp ( -i / dt'Rwit'] 

J to 



(52) 



In particular we will focus on the dependence of j(a^a) 
on the coupling strength A and the rotation velocity 6^ 
after one or many revolutions for different initial states. 
Specifically, we are interested in the effect of the time- 
dependent rotation on the mean photon number. This 
gives us an insight into the interplay of the geometric 
phase and non-equilibrium time evolution. 

In the TDL (j oo) the time evolution of the scaled 
mean photon number can be derived from the time de- 
pendent mean field equations (26l-(29). Indeed 



a^a) 



(53) 



if the system is initially (at t — to) prepared in the prod- 
uct state 



mto)) = \a{to)mto)) 



(54) 



of the field coherent state \a) and the spin coherent state 

10, 



Ht)) 



e 2 



n=0 



1 



m+j 



(55) 

\j,in). (56) 



2. Parity operator 

An additional information on dynamical behavior is 
provided by the observable, which has a meaning of dis- 
crete topological number, namely the parity of a time- 
evolved state. The parity operator II of the Dicke model 
reads 



n 



cxp 



where 



N = a^a + +j 



(57) 



(58) 



counts the number of excited quanta in the system. Ap- 
plying it to a basis state |7i)|j, m) yields 



fl|n)|j,m) =e"("+'"+^V)|j,"^), 



(59) 



since a^a|n) = n\n), J^jj,™) = m\j,m) and J^|j,m) ~ 
j{j + l)\j,m). Thus the parity operator has the eigen- 
states \n)\j,m) and the corresponding eigenvalues ±1 
(depending whether n + m + j is even or odd) . 



The time evolution of the parity operator is given by 

oo j 

mmm)-J2 E i'«".™wi'e"("+™+^"\ (6o) 

n— rn ——j 

where Kn,m{t) is defined by the expansion of a state \'ip{t)) 
in the basis \n)\j,m), i.e.. 



n— m——j 



(61) 



In the TDL {j — > oo) we use the solutions of the mean 
field equations (l26|-([29l) 



V2 V4j- 



(62) 



to calculate the time evolution of the parity for the 
dynamics starting from a coherent state \ip(to)) — 
|a(to))|C(*o))- The time evolution of the parity can be 
calculated by using the following properties of the field 
coherent states 



13 e"a. 



(a|/3) = exp 
exp(A.ata) = ^^^— ^(at)'=(a) 



k=0 

and the spin coherent states 



77 := e'-C, 



(C\v) = ^ ' 



(63) 
(64) 

(65) 



(66) 
(67) 



The time evolution of the parity in the coherent state 
\a)\() therefore reads 



(«|(C|n|C)|a) =exp(aa*(e--l)) (^^^ 



i-CC 



2j 



(68) 



3. Numerical procedure 

For a finite number of two-level atoms TV — 2j or 
pseudo-spin of finite length j = N/2 the time evolu- 
tion can not be calculated exactly. In addition, to study 
the validity of the time dependent mean field evolution 
equations we compare the mean field solution with exact 
numerical solution of the Schrodinger equation. For the 
numerical calculations we truncate the bosonic Hilbert 
space up to tim bosons but we always keep the full 
Hilbert space of the pseudo-spin. We kept um always 
high enough to assure that the error of the numerical 
data is on the level of the machine precision. For the 
numerical solution we use the Chebyshev scheme [531 to 
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calculate the time evolution. First, we apply a transfor- 
mation into a co-rotating basis |'!/')rot = R-z(^)IV')' 



(V'(t/)|aTa|V(i/)) 



ROT 



|V'(io))ROT- 

(69) 



Then we compute the time propagator of the time inde- 
pendent Hamiltonian 



H 



ROT 



(ljq + 5^)Jz + wa' a + 



A 

M 



(at + a) (j 



which is then expanded in a Chebyshev series: 

M 



-iHrotA* 



OfcTfe (hROT), 



(70) 



(71) 



fc=0 



where T'fc(hROT) are the Chebyshev polynomials of order 
k and Iirot is a rescaled Hamiltonian 

{- 2HroT — (^^Max + -^Min)! /„„\ 

ilROT - ^ • ['^) 

-C/Max — -C'Min 

Here i?Max and i?Min are the largest and the smallest 
eigenvalue of Hrqt- The expansion coefficients are 
determined by 



ak =(-«)"■ cxp ( -iAt^iEuiix + Euin) ] (2 - Sk,o) 



>/fc(Ai-(£'Max — -^Min)), 



(73) 



where ^fe^o is the Kronecker delta and Jk{x) are the Bessel 
functions of the first kind. 

In the same way we also compute the time evolution 
of the parity operator for a system with a finite number 
of two-level atoms. 



B. Stationary Dicke initial state: |acst)ICcst) 

First, we study the mean photon number if the system 
is initially (at to = 0) prepared in the product state of 
the field and spin coherent states 



with 



2A 



0, 



IV'(O)) = ® \CcJ 



,0, A < iy^wo' 



(74) 



(75) 



(76) 



Coherent state with these parameters corresponds to 
the stationary state of the time independent (unrotated) 



Dicke model (equations (26l-(29) for 5^ = 0). This can 



be seen by inverting the parametrization used previously, 
see Eq. (20) and Eq. (21). Therefore, we obtain 



91(0) 



) , A > IV^wo 



(77) 



,0, A < i^wwo 

pi(0)=0, (78) 

,,(0) ^^i^, (79) 

[O, A < ^y/UJUJo' 

P2(0)=0. (80) 

Our driving protocol is the following: we start in the ini- 
tial state IcKcstllCcBt) and then we switch on the rotation 
at to = and let the system evolve until tf — npj 27r/(50, 
where tir indicates the number of rotations. In order 
to understand the influence of the rotation we also com- 
puted the time evolution where we do not switch on the 
rotation at tg = but where we let the system evolve 
according to the usual Dicke Hamiltonian Hq. 

1. Time dependence of the mean photon number 



rot: a^=1.0, i = 10, 11^=60 
no rot : i^=0, j = 10, iij,=60 



rot : aj = 1.0, TDL 
no rot: S.=0, TDL 




Figure 3: Time dependence of the scaled mean photon num- 
ber on resonance uj = uq — 1.0 for different values of the cou- 
pling strength A = 0.6, 1.0, 1.3. The red and cyan curves are 
the time evolution obtained from the time dependent mean 
field equations and the blue and green from the numerical 
integration of the time dependent Schrodinger equation for a 
finite sized system. Here the red and blue lines correspond to 
rotationally driven system with a driving velocity of = 1.0 
and the cyan and green curves correspond to the time evolu- 
tion with no rotational driving. 



Fig. [3] compares the time evolution of the scaled mean 
photon number, {a,^a){t)/j for the rotationally driven 
Dicke model and the time evolution of the undriven Dicke 
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model. Part of the data refer to the TDL at j — )■ oo, using 
the mean field solution (red and cyan), while the other 
part corresponds to the system with a finite number of 
two level atoms (blue and green). We observe that the 
mean field solution works very well in the strong coupling 
regime and it becomes more and more accurate with in- 
creasing j, as expected. 



2. Dependence of the mean photon number on the 
atom-field coupling strength \ 

The dependence of the scaled mean photon number 
on the coupling strength A after time evolution of the 
system corresponding to a single circle, i/ = in/ 6^, is 
shown in Fig. |4]for a fixed driving velocity = 1.0 and 
on resonance oj = loq = 1.0. Again, the mean photon 

w = 1.0, aj(, = 1.0, (5^ = 1.0, tf = 2TT/S^ 
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2.0 

+ 

<0 1.5 





no rot: 5^—0, — 70 
rot: (5^ — 1.0, j=Q, n^j = 70 
rot: 5^ = 1.0, TDL 

\aj\ 5,^ = 1.0 
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A 



Figure 4: The scaled mean photon number as a function 
of the coupling strength. Comparing the mean field solution 
(j — >■ oo) with the numerical solution for a system with a 
finite number of atoms. We also compare the rotationally 
driven Dicke model with the undriven Dicke model. 



number is computed using the time dependent mean field 
equations which should be valid in the TDL (j — >■ oo, 
see orange curve) and also by solving the Schrodinger 
equation numerically (here j = 6, see green curve). We 
also compare the results with the evolution of the Dicke 
model in the absence of rotations (blue curve). 

A clear feature of the rotational driving is the appear- 
ance of a minimum in the super radiant phase (in the case 
of the Fig. |4] around A ~ 0.82). By looking into the time 
evolution of the scaled mean photon number we observe 
from Fig. [3] that the time evolution is not periodic in the 
region of the minimum which may be attributed to the 
onset of chaotic behavior noticed in Ref. [18]. The new 
minimum is an indication of the competition between two 
stationary states, corresponding to driven and unrotated 



Dicke models. 

Now let us consider an interesting limit of multiple 
rotations 0/ ^ 1. This limit is interesting because of 
a competition between dynamical and geometric phases 
|14) . In this limit it is natural to look into the time 
averaged scaled mean photon number 

((ata))T/j dt (ata)/j. (81) 

In Fig. [5] we illustrate the time averaged scaled mean 
photon number obtained for a fixed rotation velocity 8^ = 
1.0 by solving the mean field equations (orange curve) 
and by solving the Schrodinger equation for a system 
with j = 6 (green curve). We conclude therefore that 



w = 1.0, aj(, = 1.0, (5^ = 1.0, t^ = 150-27r/($^ 




Figure 5: Time averaged mean photon number as a function 
of A. The presence of a new minimum in the super radiant 
phase indicates a reentrant behavior as a consequence of a 
competition between geometric and dynamical phases |14j . 

the averaging over many rotation reveals a meta-stable 
phase in the super-radiant phase A > Ac. This phase is 
determined by a dynamical critical atom-field coupling 

strength A^'^^"''. In the case of cj = = 1-0 and 5^ ~ 1.0 
it is given by Ai'^^"^ « 0.82. 

3. Dependence of the mean photon number on the rotation 
velocity 

The dependence of the scaled mean photon number 
on the driving velocity after evolving the system for one 
circle tf = 2i:/5^, with different atom-field coupling pa- 
rameters A is plotted in Fig. |6] We observe that if the 
critical paraboloid is not encircled A < Ac = ^wwo /2, 
the mean photon number is zero. But if we encircle the 
critical paraboloid the mean photon number is non-zero 
and there should be a critical driving velocity. Similar 
to the case of the stationary state of the rotated Dicke 
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= 1.0, 1^0 = 1.0, t; = 2-nl5^ 



_ rot: j =6, Km =70, A =0.3 

rot: TDL, A =0.3 

— no rot: J=6, n^i = 10, A =0.3 
- - - no rot : TDL, A =0.3 

rot: j=6, rij,=70, A = 1.0 


rot: TDL, A =1.0 

no rot : j =6, = 70, A = 1.0 

no rot : TDL, A =1.0 
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Figure 6: The driving velocity dependence of the scaled mean 
photon number after a single circle. 



model, where the critical atom-field coupling is given by 
^yLjj{(jJo + 5,f,) /2 which defines a critical driving velocity 

^lp°c^ (see the black curve in Fig. 6 \ac2 1^ / j)- Note how- 
ever that in the case of one rotation we can not see an 
indication of presence of the critical driving velocity in 
the dynamics. 

In Fig. [7] we also plot the time averaged scaled mean 
photon number which clearly illustrate the presence of a 
dynamical or rotational critical driving velocity S^l^^'^^ in 
the TDL. 

aj = 1.0, W(, = 1.0, A=1.0, t^ = 150-27r/J^ 



0.0 






































































rot : j = 6, ri^^ = 70 

- rot : TDL 

no rot: j =6, nj^^=7(i 

— - no rot : TDL 















2.0 



Figure 7: Time averaged scaled mean photon number as 
a function of the rotation velocity. We compare {{a^a))T/j 
for the rotationally driven Dicke model with the one for the 
undriven Dicke model. In both cases we kept tf = 15O-27r/(50. 



4- Non-Equilibrium Phase Diagram 

Summarizing our studies, in Fig.[8]we plot the time av- 
eraged scaled mean photon number as a function of the 
driving velocity and the atom-field coupling strength ob- 
tained from the mean field equations (26)- (29). One can 
clearly distinguish two distinct "meta stable" phases in 
the super-radiant phase. We therefore suggest this plot 
as a non-equilibrium phase diagram of the rotationally 
driven Dicke model. We associate this reentrant behav- 
ior with a competition between geometric and dynamical 
phases. This mechanism is described in Ref. [H]. 

The line in the super-radiant phase (the red dashed 
curve in Fig. |8]) which describes a dynamical critical line. 



was fitted with the function aI'^^"'' = 0.5 



0.327 It 



defines a dynamical critical atom-filed coupling strength 
^(dyn) reciprocally a dynamical critical driving veloc- 
ity (J^^^"^ This should not be confused with the criti- 
cal coupling Ac'°*^ — y^uj{ujo + 6^) '/2, that emerges if we 
start the evolution of a system from the stationary state 



|ac2)|Cc2) (see next paragraph). 




= 0)0 = 1, tj: = 150-27r/(5^, in TDL 



A 0.7 



A = 0.5 




Figure 8: Dependence of the time averaged scaled mean 
photon number on velocity of rotation and the atom-field cou- 
pling strength. It is obtained from the mean field equations 
for ij — Ldo ~ 1.0 and averaged over a time tf = 150 ■ 27r/(5^. 
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5. Parity Yl 

In Fig. |9]we show the time evolution of the parity oper- 
ator starting in the coherent state |ast)|Cst) for a system 
with a finite number of two-level atoms. The parity is 
constant in time and once we enter a super-radiant phase 
A > Ac = ^wwq /2 it starts to drop from one to zero. The 
inset plot shows the A dependence of the parity and its 
tendency to a step function with growing j . 



1.5 



1.0 



a; — 1.0, oJq — 1.0, (5^ — 1.0, 



= 100 



0.5 - 



<U 0.0 



-0.5 - 



-1.0 - 



-1.5 



J =3 

J =10 



0.4 0.5 0.6 0.7 

A 



j = 10, A =0.40 

j = 10, A=0..il 

i = 10, A =0.61 



27r 



47r Gtt 
t 



IOtt 



Figure 9: Time evolution of the parity operator when we start 
the dynamics in the coherent state |ast)|Cst) for a system with 
a finite number of atoms. The inset plot shows the depen- 
dence on the coupling constant A for dilTerent values of j. 



The time evolution of the parity (11) and its depen- 
dence on the atom-field coupling A calculated from the 
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mean field equations, i.e., by Eq. (68), is shown in Fig 
The parity is not a constant function of time in the previ 
ously observed intermediate regime of (a^a)/j. The time 
averaged parity as well reveals the observed dynamical 
critical coupling A^'^'^"^ ~ 0.82. Nevertheless, for j — oo 
the parity becomes constant in time and goes to in the 
entire super-radiant phase. 



Further, in Fig. 11 we also illustrate the scaled parity, 
namely the parity where all the phase space coordinates 
(ft, Pi) are rescaled by : (ft, Pi) ^ {(lil \/J -.Vil \/J)- 
This implies 



(a|(c|nsclC)l«) = 



9? 



2\ 2i 

-vV 



2f 



(82) 



aj = 1.0, W(, = 1.0, Sj, = 1.0 
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Figure 10: Upper panel: time evolution of the parity oper- 
ator if we start the dynamics in the coherent state |ast)!Cst)- 
Lower left panel: the dependence on the atom-field coupling 
constant A of the parity after 20 circle. Lower right panel: the 
dependence of the time averaged parity on the coupling con- 
stant A. The parity is obtained by solving the time dependent 
mean field equations. 
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A=0.40, j=200 
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Figure 11: Upper panel: time evolution of the scaled par- 
ity operator if we start the dynamics in the coherent state 
|Q!st)|Cst)- Lower left panel: the dependence of the scaled par- 
ity on the coupling constant A after one circle. Lower right 
panel: the dependence of the time averaged scaled parity on 
the coupling constant A. The scaled parity was computed by 
solving the time dependent mean field equations. 



The scaled parity is naturally independent of the num- 
ber of spins j but it is no longer constant in time for 
couplings between 0.5 and 1. The time averaged scaled 
parity shows also an intermediate regime similar to the 
intermediate regime in the time averaged scaled mean 
photon number. 



C. Stationary circle initial state: |Qc2)|Cc2) 

Here we study the mean photon number if the system 
is initially (at tg = 0) prepared in the product state of 
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the field and spin coherent states 

|V(0)) = |«e.)®|Cc.), 

with 



etc 



(83) 

(84) 
(85) 



where = uj{ujq +(50). This choice of the parameters a 
and C of the coherent states corresponds to the station- 
ary "circle" solutions of the time dependent mean field 
equations (26l-(29). In the same way as in the previous 




section this can be obtained using Eq. ( 20 1 which gives 




(86) 
(87) 



A < 



P2(0) =0. 



(89) 



The driving protocol is the following: we start with the 



system prepared in the state 



>ICc 



Then we switch 



on the rotation at = and let the system evolve ac- 
cording to HrdI^) until = riR 27r/(S0. Here indicates 
the number of rotations in the space of parameters (see 
Fig. [ij. Further, we consider the time evolution where 
we do not switch on the rotation at io = but where we 
evolve the system according to (undriven case) in or- 
der to understand the influence of the rotational driving. 



1. Mean photon number for the stationary coherent state 

From the fixed point solution of the equations of mo- 
tions (26)- (29 1 we observe that if we start the evolu- 



tion from the state ( 83 1 , the dynamics of the system re- 
mains bounded to the circle and after a closed evolution 
tf = 2it/5^ the mean photon number in the thermody- 
namic limit is simply given by 



1 



a^a) 



1 f 2A) 

2 \ LO I 



(CcaKacala'^alacJICc 



1 

2j 



' n 







ifA>i 
if A < 



(90) 



We also observe now that the critical coupling, which 
marks the quantum phase transition, is shifted by 
the amount of the 



|\/w(wo + <50)'. In Fig. 
der to illustrate the shi: 
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(rot) 
t /\c — 

(left) we plot {a,^k)c2/j in or- 



rota tion velocity as A 



by the rotation. Further, in Fig. [12] (right) we draw 
the velocity dependence of {si^d^c^/j for circles C2 inside 
the critical paraboloid and compare it with the situa- 
tion corresponding to circles Ci that encircle the critical 
paraboloid. One can see that the critical driving veloc- 
ity is defined by 5^^°^^ = ujq. Above this criti- 
cal velocity the mean photon number is zero. This al- 
lows one to probe the equilibrium quantum critical point 
Ac = ^wcjQ /2 in an indirect way by encircling critical 
paraboloid in the parameter space and measuring the 
dependence of mean photon number as a function of a 
driving velocity. 
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Figure 12: Plot of the mean photon number (a^a)c2/j as a 
function of the coupling strength A on the left panel and as a 
function of the rotation velocity 5^ on the right panel. 



2. Deviation from stationary coherent state 

The natural question which immediately arises is what 
happens if the initial state is not exactly given by the 
fixed point conditions, namely if the initial state is not a 
coherent state 



\a{t)) = e 

m) 



a" 



n=0 
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with a(0) = and C(0) = Cca- In Fig- (13) we plot the 
scaled mean photon number obtained from solving the 
time dependent mean field equations but this time with 
the initial condition that slightly differ from the fixed 
point. 



3. Time dependence of the mean photon number 

In Fig.[T4]wc illustrate the time evolution of the scaled 
mean photon number. In the TDL (j — > oo) we used the 



t of the critical coupling caused time dependent mean field equations ( 26 )-( 29 ) in order to 
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w=aj„=l, S^ = l, j=lQQ u=w„ = l, A=l, j=lQO 




Figure 13: The mean photon number {a} a)c2l j as a function 
of the coupling strength A on the left and as a function of the 
rotation velocity 5^ on the right starting on the fixed point 
(qi,c2 ,Pi,c2 5 92,C2 ,P2,C2 ) (blue dots) and with small deviation 
from this fixed point (red + and green x). 



tj = 1.0, wo = 1.0, A = 1.0 

2.6 I I I I I 

rot ; 5^ = 1.0, j =6; n^j = 70 — no rot : 5^ =0, j =6, %j — 70 

rot: 4^ = 1.0, ;=10, nj, = 70 -- no rot : (S,j=0, j=10, Hi, = 70 

2-4- rot: 5^ = 1.0, i = 12, nj, = 70 -- no rot : (S,, =0, j = 12, rij, = 70 

rot: (S^ = 1.0, TDL no rot : (S^ = 0, TDL 




compute {a^a){t)/j (orange curve) while for a finite num- 
ber of two- level systems we used the Chebyshev technique 
to calculate the time evolution. We compare the time 
evolution of the scaled mean photon number of the rota- 
tionally driven Dicke model (solid lines in Fig. 14) with 



the time evolution of the scaled mean photon number 
of the Dicke model with no rotation of the pseudo-spin 
vector around the z axis (dashed lines in Fig. 



14) 



As expected, if we start in the coherent state |ckc2)ICc2) 
the time evolution of {a}a}/j is constant in the TDL 
(solid orange line in Fig. 14). However, for finite j one 



can see that the scaled mean photon number is not time 
independent, but increasing j it tends to a constant value 
(the amplitude of the oscillation decreases). In the case 
where we let the system evolve without rotational driv- 
ing, the scaled mean photon number shows pronounced 
oscillatory behavior compared to the rotational driving 
which is not completely 27r periodic. 



4- Dependence of the mean photon number on the 
atom-field coupling strength A 

In this subsection we illustrate the dependence of the 
scaled mean photon numbers on the atom-filed coupling 
strength, (a^a)(A)/j. The system is initially prepared in 
the product coherent state |ac2)|Cc2)- First we consider 
only one revolution in parameter space, tf — l-n jb^. For 
rotational driving the system remains in the same state 
and there is no significant deviation from the static be- 
havior of the scaled mean photon number. However, if 
the system is evolved according to the usual Dicke model 
Hd a peak is observed at . For A < Ai"^"*^ there is a 
large discrepancy between the finite size calculation and 
the mean field solution. 

In Fig. [16] we show the time averaged s cale d mean pho- 
ton number, ((a^a))7-/j as defined in Eq. 81 for the fixed 
rotation velocity 5^ 



Figure 14: Time dependence of the scaled mean photon 
number, (a^a)/j, for different values of j, a coupling A — 1.0 
and for to = ijjq = 1.0. The solid lines represent the time 
evolution with rotational driving and the dashed lines give the 
time evolution of the Dicke model without rotational driving. 
The evolution time tf m. both cases is the same. 



rot: (S,j = 1.0, j=6, n^i=10 
rot: 5^ = 1.0, ] = 12, nn = 70 
rot: <S^ = 1.0, TDL 
no rot: 5,(,=0, ^'=6, /i,^j = 70 



no rot: 5,(,=0, ^' = 12, /i-^j = 70 
no rot: 5^=0, TDL 




Figure 15: The scaled mean photon number as a function 
of the coupling strength A. We compare {a^a)/j for the rota- 
tionally driven Dicke model with the one for undriven Dicke 
model {tf = 2ti/S^ was kept the same in both cases). 



level systems N — 2j the time propagation was calcu- 
lated using the Chebyshev scheme [24] and in the TDL 
we solved the time-dependent mean field equations. 

For the initial state |q;c2)ICc2) which corresponds to 
the stationary state of Hrd we effectively observe no dif- 
ference between \ac2\ /j and the time evolved compu- 



1.0. For a finite number of two- tations. However, when the system is evolved accord- 
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w = 1.0, a;o=1.0, (5^ = 1.0, ty = 150-27r/<5( 



• — ■ rot: i^ = 1.0, ;=6, iij, = 70 
4.0 - rot: (S,j = 1.0, TDL 



no rot : 5^ = 0, j =6, = 70 




Figure 16: The time averaged scaled mean piioton number as 
a function of the couphng strength A. We compare ({a^a))T/j 
for the rotationally driven Dicke model and the unrotated 
Dicke model {tf was kept the same in both cases). 



uj = l.Q, a;„ = 1.0, tj = 2TT/S, 
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Figure 17; The scaled mean photon number as a function 
of the rotation velocity 5^. We compare (a^a)/j for the rota- 
tionally driven Dicke model and for the undriven Dicke model 
with tf =2n/5^ in both cases. 



ing to Hd we again observe that below A^"^"'^ the finite 
size results do not coincide in both cases. This can be 
corrected if we introduce small fluctuations to the semi- 
classical description by setting the initial conditions of 
the mean field equations ([26|-(29) not strictly to zero 
but a(0) = C(0) = a*(0) = C*(0) = 10"^ (brown line). 



5. Dependence of the mean photon number on the rotation 
velocity 5^ 

In this subsection wc illustrate the dependence of the 
scaled mean photon number on the rotation velocity (5^ 
after one revolution in parameter space, tf = 2tt/5^^ 
which is shown in Fig. 17 We observe that if the critical 
paraboloid is encircled, the scaled mean photon number 
is different from zero and the critical driving velocity is 

given by (J^"*'' = ^ uJo- This basically follows from 

the rotation velocity shift in the rotated critical coupling 

,(rot) 



Ac'"''' — Y^aJ(c[Jo~+~(50)~'/2. This allows one to probe the 
equilibrium quantum critical point Ac = ^JljlJq /2 from a 
distance in parameter space. 

The velocity dependence of the time averaged scaled 



mean photon number (Fig. 18) clearly reveals the crit- 
ical driving velocity i^^"*'*. However, if we suddenly, at 
^0 = 0, switch on the undriven evolution (with Hq) the 

mean field solution is zero for (5^ > ^0°*'' while the nu- 
merical calculation for a finite number of atoms is non 
zero in this region. Again, this can be corrected if we al- 
low small fluctuations in the semi-classical description by 
setting the initial conditions of the mean field equations 
not strictly to zero: a(0) = C(0) = a*(0) = C*(0) = 10"^ 



in the corresponding range of the parameters. 

[j = 1.0, ^0 = 1.0, A=1.0, ty = 150 -2^/(5^ 
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Figure 18: The time averaged scaled mean photon number as 
a function of the rotation velocity (5^. We compare ({a^a))T/j 
for the rotationally driven Dicke model with the results for the 
undriven Dicke model averaged over tf — 150 • 2t7/S^ in both 
cases. 



6. Non-Equilibrium Phase Diagram 

We summarize the dependence of the time averaged 
scaled mean photon number on the atom field-coupling 
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strength A and the rotation velocity S^f, in Fig. 19 Here 
we see that the rotational driving shifts the super-radiant 
phase transition by the amount given by the rotation 
velocity. The critical line (solid red line in Fig. 
defined by A^™" = ^ujicjo + 5^y or by S^lf = ^ 



19|) is 

-Wo- 



A 0.7 



A = 0.5 




Figure 19: The time averaged scaled mean photon number as 
a function of the rotation velocity and the atom-field coupling 
strength. Results are obtained from the mean field equations 
for Lu = coo = 1-0 and averaged over a time tf = 150 • 2-k/S^. 



7. Parity H 

The time evolution of the parity operator 
is calculated for the stationary circle initial state 

>^C2)ICc2) 
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for a system 
Note that the 



lo^ca^Kca)- Results are shown in Fig. 
with a finite number of two-level atoms, 
parity is constant in time and jumps from one to zero at 
the critical atom- field coupling strength Ai™*-*. 

In the TDL [j oo) we use the solutions of the mean 
field equations to calculate the time evolution of the par- 
ity, Fig. [21] We also illustrate the time evolution of the 
scaled parity (where as previously we rescaled all the 
phase space coordinates by vT)- 
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Figure 20: Time evolution of the parity operator starting 
from the coherent state |Qc2)ICc2) for ^ system of finite number 
of two-level atoms. The inset plot shows the dependence on 
the coupling constant A for different values of j. 
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Figure 21: Upper panel: time evolution of the parity (11) 
starting from the coherent state |Q!c2)ICc2)- Lower left panel: 
dependence of (11) on the coupling constant A for different 
values of j after 50 circles. Lower right panel: time averaged 
parity as a function of the coupling A. The time evolution is 
calculated from the time dependent mean field equations. 



D. 



Fock initial state: \n = 0)\j,m — —j) 



Here we consider the evolution of the mean photon 
number if we initially prepare the systems in the Fock 
state with no photons n = while the spin state has the 
lowest weight, m = —j, 



(93) 
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Figure 22: Upper panel: time evolution of the scaled parity 
operator starting from the coherent state |Q!c2)ICc2)- Lower 
left panel: the dependence on the coupling constant A of the 
scaled parity. Lower right panel: the time averaged scaled 
parity as a function of A. 



Figure 23: Time dependence of the scaled mean photon 
number on a resonance lj = ljq = 1.0 for different values of 
the coupling strength A. The system is initially prepared in 
the state |0)| j, — j). 



We switch on the rotation at = and let the system 
evolve until tf = riR I'k jb^^ where we measure the mean 
photon number {a}a){tf)/j. We also compare it with 
the evolution of the un-rotated system. In terms of the 
coherent states this initial state corresponds to a = 
and C = 0- The mean photon number is calculated in the 



TDL using the mean field equations (26)-(29l with the 
initial conditions a(0) = a*{0) = C(0) = C*(0) = 0. 



1. Time evolution of the mean photon number 



In Fig. 23 we show the time evolution of the scaled 
mean photon number, (a^a)/j. It compares the scaled 
mean photon calculated for a system of finite size (j = 
6) by the Chebyshev scheme with and without rotations 
(red and blue curve respectively) with results for the TDL 
(j — >■ oo) with and without rotations (orange and green 
curve respectively). 

In the static case the mean photon number is explic- 
itly zero for this initial state, however for the driven situa- 
tion we observe a critical coupling above which the scaled 
mean photon number becomes macroscopic. We observe 
that for the initial Fock state, \0)\j, — j), the scaled mean 
photon number in the rotated case (we switch on the ro- 
tation at ^0 = 0) is still microscopic for A = 0.5 whereas if 
no rotation is switched on the scaled mean photon num- 
ber becomes macroscopic at Ac = ^wwq '/ 2 = 0.5. There- 
fore, the critical coupling Ac is shifted by the amount 
given by the applied rotation velocity (5^. Thus for the 

driven case Ai'°'^ = y^uj{LUQ + S^j,) /2. 

Note that in the TDL the mean field equations do not 
reproduce the finite size calculations of the scaled mean 



photon number. 



2. Dependence of the mean photon number on the 
atom-field coupling strength A 

Accordingly, we study the dependence of the scaled 
mean photon number on the atom-field coupling strength 
A. First only for one closed circle in parameter space 
(/)/ = 27r, i.e., tf = 2^/6^. The scaled mean photon num- 
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Figure 24: Scaled mean photon number dependence on the 
atom- field coupling strength on resonance ui = = l.Q for a 
rotation velocity 5^ = 1.0. The system at to — is prepared 
in the state \n = 0)\j, m = 
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ber becomes macroscopic for a coupling strength greater 
than some critical coupling. In the unrotated case this 
is given by Ac = ^cjwq /2 while for a rotational driving 

it is given by A^™'^ = ^Juj(wq + 5^) '/2. In the TDL the 
mean field solutions give zero. 

By looking at the time evolution of the scaled mean 
photon number one can observe rather irregular structure 



of the oscillations. Motivated by this in Fig. (25) we also 
plot the time averaged scaled mean photon number as a 
fimction of A for multiple rotations t/ = 150 -271/(50. This 

u) = l.O, u)Q = l.O, (5^ = 1.0, tj, = 15O-27r/'50 
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no rot: 5^=0, j=6, n^y = 70 
no rot: i^=0, j = 10, i!j, = 100 
no rot : 6^=0, TDL 
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Figure 25; Time averaged scaled mean photon number on 
resonance uj — ljo ~ ^-0 for a rotation velocity 5^ — 1.0. The 
system at to = is prepared in the state \n = 0)\j, —j) and 
evolved until tf = 150 ■ 2tt/S^. The mean field solutions in 
the TDL remain strictly zero. 



plot shows that if we prepare the initial state to be a Fock 
state \n — 0)|j, to = —j) and evolve the system in time 
we observe two distinct phases. Below a critical coupling 
strength the time averaged scaled mean photon number 
is zero (normal phase) while above a critical coupling the 
time averaged scaled mean photon number becomes non- 
zero (super-radiant phase). Rotational driving shifts this 

critical coupling strength as Xc°^^ = ^Ju}(ojq + S^) /2. 
We therefore conclude here that the evolution starting 
from this state is similar to the one corresponding to the 



We observe that if we encircle the critical paraboloid (C2) 
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Figure 26: Dependence of {a^a)/j on rotation velocity 5^ 
after one circular evolution in the parameter space for a fixed 
atom-field coupling A = 1.0. The mean field solution does not 
reproduce the finite size calculations it is strictly zero. 



the scaled mean photon number shows velocity depen- 
dence. Fig. [27] shows the time averaged scaled mean 
photon number {{s) a))T / j as a function of the rota- 
tion velocity 5^. It was averaged over a time interval 



of tf = 150 -271/(50. 
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3. Dependence of the mean photon number on the rotation 
velocity 5^ 

Here we look into the dependence of the scaled mean 
photon number on the rotation velocity (5^, {d? d){5^) / j 
when the system is initially prepared in the Fock state 
IV'(O)) — |n = 0)|j, TO = — j). Again first, we consider 
only one rotation in the parameter space, tf = 21x18^ 
and compare it to the undriven evolution (see Fig. |26[). 



Figure 27: Dependence of the time averaged scaled mean 
photon number on the driving velocity at resonance tti = wo = 
1.0 and for A = 1.0. At = the system is prepared in the 
state \n = 0)|j, m = — j) and evolved up to tj = 150 ■ 2-k/5^. 
The mean field solutions always stays zero. 

It is clear that there is a critical driving velocity for 
rotational evolution given by 5^^°^^ = ^ ujq allow- 
ing to probe two different non-equilibrium phases of the 
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Dicke model. On the other hand, in the case of the time 
evolution of the undriven Dicke Hamiltonian the time av- 
eraged scaled mean photon number becomes constant (if 
we encircle the critical paraboloid). This shows that the 
presence of the non-equilibrium transition in the Dicke 
model can depend on the initial state. 



J^. Parity tl 

Here we study the time evolution of the parity operator 
if the system is initially prepared in the Fock state. In 
Fig. [28] we show the time evolution of the parity operator 
(n) for different coupling strength at fixed finite spin 
j — 10. In this case the parity is always constant in time 
with no sign of change to zero at a critical coupling. The 
system always remains in the state with the same parity 
+1. 
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Figure 28: Time evolution of the parity operator when the 
system is driven from the initial Fock state \ n = 0)|j, m = — j) 
for a pseudo spin of j = 10. The inset shows the dependence 
on the atom-file coupling strength A for different pseudo-spin 
length j. 



In the TDL (j — > oo) one can use the solutions of the 
mean field equations. In this case the parity is also un- 
changed (see Fig. 29 ) which can be traced back to the fact 
that the initial conditions a(0) = a*(0) = C(0) = C*(0) = 
allow no quantum fluctuations and therefore do not 
evolve in time. Also in the case of the scaled phase space 
coordinates (91,7)1,(72,^2) l/vT('?i,Pi, 92,^2) the par- 
ity is always +1. 
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Figure 29: Upper panel: time evolution of the parity oper- 
ator for the Fock initial state \n = 0)| j, m — — j). Lower left 
panel: dependence of the parity on the coupling constant A 
after 20 circles. Lower right panel: dependence of the time 
averaged parity ((n))T on the atom- field coupling A. 



by the Fock state with no photons n = and lowest 
weight spin state m = —j, |V'(0)) — \n — 0)|j,m = —j). 
This correspond to zero parameters of the coherent states 
a = and C = 0. This state does not evolve in time in the 
semi-classical limit. As shown in the previous paragraph 
the mean field equations did not reproduce the results 
calculated numerically for a system with a finite number 
of two-level atoms. 

Motivated by this we study the effect of small pertur- 
bation of the initial conditions. We consider a situation 
when the parameters a and ^ are not strictly zero but 
small and parametrized as 



a(0) = 10" 



C(0) = 10- 



(94) 



where e is large. Consider the following initial coherent 
state 

mo = 0)) = |a(0) = 10-^) ® |C(0) = 10-'), (95) 

which differs slightly from the Fock state with no 
photons and lowest weight state: for example, for 
j ~ 10 and hm — 100 we find for the overlap 
{a = 10-^C = 10"^ |n = 0,™ = -j) = 0.99999. 

We start the non-equilibrium evolution by switching 
on the rotation at tg = and evolving the system until 
tf = '^R27r/(50. Again we will compare the rotationally 
driven evolution and the free evolution. 



E. Nearly Fock initial state: 

|a(0) = 10-=)lC(0) = 10-^) 



Time evolution of the mean photon number 



In the previous subsection we looked into the rotational 
dynamics of the Dicke model with the initial state given 



The time evolution of the scal ed mean photon number, 
{a,'^k)/j, is depicted in Fig. (30). It compares the scaled 
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mean photon number obtained for a finite system using 
the numerical Chebyshev scheme (blue and red curve) 
and the results of the TDL computed from the mean 
field equations (26)-([29| (green and orange curve). We 
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Figure 30: Time dependence of the scaled mean photon 
number on resonance u — uuo ~ 1-0 and A = 1.0 for two 
different rotations velocities: one is slightly below the criti- 
cal driving velocity S^^c = 3.0 while the second one is just 
above it. The system is initially prepared in the coherent 
state |a(0) = 10"^)|C(0) = 10"^), with e = 3 and e = 6 for 
rotational driving. 



observe that the initial state |a(0) = 10"')|C(0) = 10"') 
for different e does generate the same dynamics as the 
initial Fock state \n = 0)|j, m = — j) for a finite number 
of two-level atoms. However, in the TDL the mean field 
equations show macroscopic excitations if we encircle the 
criticality (here A = 1.0 > \c — 0.5) and drive the system 

slower than the critical driving velocity 8^^c = wo. 

We note that the "soliton train" appearing in this figure 
(green and yellow lines) is somewhat reminiscent to the 
findings of Ref. [15] , HZ] • Note however that the model 
we study here is non-integrable. 



In Fig. (31 1 we show the time evolution of the scaled 



mean photon number with no rotational driving, i.e., the 
time evolution is governed by the usual Dicke Hamilto- 
nian Hd. Again, the system is initially prepared in the 
coherent state |a(0) = 10-'^)|C(0) = 10"'). 

Summarizing these results we observe that the initial 
state |a(0) = 10~')|C(0) = 10~') for different e and finite 
J reproduces the same dynamics as for the initial Fock 
state \n — 0)|j, m = —§)■ In the TDL the mean field 
equations show a macroscopic excitations of the mean 
photon number in the super-radiant phase (A > Ac) if we 
use the initial conditions a(0) = a*(0) = C(0) = C*(0) — 
lO"*^. This means that small fluctuations present in the 
initial state make the quasi-classical description possible. 
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no rot: j=10, 71^ = 100, e=6 



no rot: TDL, e=3 
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Figure 31: Time dependence of the scaled mean pho- 
ton number on resonance = ijJq = \X) for two differ- 
ent coupling strength A: one slightly below the critical 
driving velocity Ac = = 0.5 and one just above 

it. The system is initially prepared in the coherent state 
|a(0) = 10"')|C(0) = 10"'), with e = 3 and e = 6 and evolved 
with undriven Hamiltonian. 



2. Dependence of the mean photon number on the 
atom-field coupling strength A 



In Fig. 32 



we plot {{a}a))T/j as a function A for the 
driven model (blue, red, green) and for the undriven one 
(cyan, orange, lime). We again observe a critical coupling 
strength Ac at which the number of excited photons be- 
comes macroscopic. In the case of rotational driving the 
critical coupling is shifted by the amount of the applied 
rotation velocity S^. The A-dependence of {{a^k))T/j for 
the initial Fock state of the Dicke model with a finite 
number of atoms looks similar to the one obtained in the 
TDL with a(0) = a*(0) = C(0) = C*(0) = 10-^ How- 
ever, the main differences with a driven model are (i) 
appearance of dips and (ii) difference in the scaling form 



(rot) 



of the curve when A — >■ A 



3. Dependence of the mean photon number on the rotation 
velocity 5^ 

The dependence of {{a^a))T/j on for a driven (blue, 
red, green) and undriven (cyan, orange, lime) systems is 
shown in Fig. p3| ). When the critical paraboloid (path 
C2) is encircled, there is a critical driving velocity given 

by (5^™*^ — ^ 1^0- This is in contrast to the TDL with 

the initial conditions a(0) = a*(0) = C(0) = C*(0) = 
when the mean photon number was always zero. In the 
case of undriven model the time averaged scaled mean 
photon number fiuctuates around a constant value. 
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Figure 32: Time-averaged mean photon number ({a^a))T/,7 
as a function of the coupling strength. In the TDL (j — )■ cxd) 
obtained by solving the mean field equations (26l-(29l with 
the initial conditions q(0) = q*(0) = C(0) = ('(0) = 10"'. 



Figure 34: Upper panel: time evolution of the parity opera- 
tor for the driven system starting from the nearly-Fock state. 
Lower left panel: the parity after 50 circles as a function of 
the coupling. Lower right plot: the same as left panel but 
averaged over 150 circles. 
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Figure 33: Time-averaged mean photon number {(a^a))T/j 
as a function of the driving velocity. The TDL (j — )■ oo) is 
obtained by solving the mean field equations (|26[)-(|29[) with 
the initial conditions q(0) = a*(0) = C(0) = C*(OT= 10"'. 



4- Parity II 



a smooth jumps from +1 to 0. We also show the par- 
ity calculated using the rescaled phase space coordinates 
(see Fig. 35 1. 
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Using the initial state |a(0) = 10-')|C(0) = 10"') with 
e > 2 for a system of a finite number of atoms we find 
the same behavior, (11) = +1, for any time and coupling 
strength as for the initial Fock state \n = 0)\j,in — — j). 
On the other hand, in the TDL {j oo) using the so- 
lutions of the mean field equations (26)-([29| we compute 



the time evolution of the parity and observe (Fig. 34) 
that for a coupling A > Ac the time evolution shows 



Figure 35: Upper panel: time evolution of the scaled parity 
operator when we start the rotational dynamics with the sys- 
tem prepared in the nearly-Fock state. Lower left panel: the 
scaled parity as a function of A after 50 circles. Lower right 
panel: the same as the left one but averaged over 150 circles. 
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F. Ground state of the un-rotated Dicke model as 
initial state: |GS) 

Here we consider the evolution of a system initially 
prepared in the ground state 

|V(to = 0))-|GS). (96) 

First, we numerically compute the ground state |GS) 
of the Dicke Hamiltonian Hd by truncating the bosonic 
Hilbert space up to um while maintaining the full Hilbert 
space of the pseudo-spin j. Then, at to = 0, we start to 
evolve the system according to the rotated Dicke model 
HRD(i)- Next we measure the scaled mean photon num- 
ber (a^a) as 



(at a) 



= (GS|e™"°^*^ atae-*'^«°^*^|GS), 



(97) 



where the final time is, as before, given by i/ = tt-r 2tt/6^. 
We also consider the undriven evolution governed by the 
evolution of the usual Dicke model Hp. 



1. Time evolution of the mean photon number 

The time dependence of the scaled mean photon num- 
ber is shown in Fig. [36] Here we compare the results for 
the initial ground state |GS) (red curves) with the results 
for the initial coherent state with the parameters ac^^ and 
Ccrt (blue and green curves). The dashed lines show the 
situation for the undriven evolution. 
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Figure 36: Time dependence of the scaled mean photon 
number on resonance uj = luq — 1.0 for different values of the 
coupling strength A. Here we compare the results for different 
initial states: |GS) and |Qc3t)|Ccst)- The ground state |GS) is 
obtained numerically and we choose um ~ 100. 



ber. For a finite number of atoms the time evolution of 
the scaled mean photon number for the ground state | GS) 
is close to the one for the coherent state lac^t) ICc^)- Foi" 
large atom-field couplings they become identical. Thus, 
the coherent state |acgt)lCcat) ^ good approximation 
for the ground state of the system even in the driven 
case. This phenomenon is noticed already in |22| . In 
the undriven case the time evolution of the mean photon 
number for the ground state is constant and corresponds 
to the value found in the TDL 



(GS|aU|GS) = 



1 (2X^ 

2 \ cj J 



1 



\ 4A2 ) 
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2. Dependence of the mean photon number on the 
atom-field coupling strength A 

Here we study the dependence of the scaled mean 
photon number on the atom-field coupling strength A. 
First, we consider one revolution in parameter space, 
tf — 27r/(50. This case is shown for a fixed rotation ve- 
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locity 6^ — 1.0 at resonance oj — (jJq = 1-0, in Fig 
We also show the dependence of the scaled mean photon 
number on A for the undriven Dicke model (green curve) , 
where the final time was kept the same tf = 2Tr/S^. 

aj = 1.0, Wo = 1.0, <$^ = 1.0, tf = 2TT/S^ 
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Figure 37: The scaled mean photon number as a function of 
the atom-field coupling strength: here we compare the mean 
field solution (j — ^ oo) for the initial coherent state |Qcat)ICcat) 
(yellow solid line) with the time evolution calculated numer- 
ically for a system with a finite number of atoms {j = 6) for 
the initial state chosen as ground state |GS) (blue line) and 
the coherent state (red fine). We also compare it with the 
undriven case (green fine) and the static (dashed lime). 



As we have already seen the time dependence of the 
For the rotationally driven case (solid lines) we observe scaled mean photon number calculated for the initial 
an oscillatory behavior of the scaled mean photon num- ground state |GS) and for the initial coherent state 
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Ittc^JICcst) give qualitatively the same results. For the 
ground initial state we observe a minimum in the super- 
radiant phase (in Fig. 37 at A « 0.82). This indicates 
a competition between the two stationary states or the 
ground states of the rotationally driven and undriven 
Dicke Hamiltonians. One further interesting point is the 
difference in scaling behavior close to the critical point 
Ac = ^wwq '/2: This behavior is different for the ini- 
tial ground and coherent states. For the ground state 
it follows the same scaling behavior as in the static and 
unrotated case. 

So far we considered the mean photon number after 
one circle = 27r, but the situation with many rota- 
tions 0/ — ?> oo represent an interesting limit. In this case 
the interplay of dynamical and geometrical phase may be 
a source some new effects p3]- Considering the case of 



many rotations in Fig. 36 it is natural to take the time 
average of the scaled mean photon number as already 
defined in Eq. (81 ). The time averaged scaled mean pho- 

Here we 
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ton number {{aJa))T/j is depicted in Fig. 
averaged over a time interval of — 150 • 2Tr/S^ for the 
resonance conditions uj — ujq = 1.0 and kept the rotation 
velocity to be — 1.0. 

tj = 1.0, ^0 = 1.0, (5^ = 1.0, t^ = 150-27r/<5^ 



GS), rot: i^ = 1.0, j=6, n„=70 

™*- *»=1-''' J=S' "11='''' 
|a^C,,>, r"t: *8=l-0, TDL 
|GS), no rot: S^=0, j=6, n^=7Q 




Figure 38: The time averaged scaled mean photon number 
{{a^a))T/j at resonance uj — ujo = 1.0 as a function of the 
coupling A. The system is initially prepared in the ground 
state |GS) which is calculated numerically and then evolved 
until tf = 150 • 2n/5^ (blue curve). The result is compared 
with the one obtained for the initial coherent state (red and 
yellow curve) and with the one corresponding to undriven 
system (green line). 



We qualitatively observe the same behavior for the 
time averaged scaled mean photon number for the two 
different initial states, namely for the ground state |GS) 
and for the coherent state |acrt)ICcst)- Moreover, the time 
averaged scaled mean photon number shows a different 
scaling for the two initial states when approaching the 
equilibrium critical coupling Ac = ^wojq '/2. 



Dependence of the mean photon number on the rotation 
velocity 5^ 



In this subsection we illustrate the dependence of the 
scaled mean photon number on the rotation velocity S^. 
First, in Fig. 39 we consider one revolution, tf = 2Tr/5^ 
for different atom-field couplings and for ui = = 1.0. 
We plot the velocity dependence of the scaled mean pho- 
ton number for rotational driving once for the initial 
ground state |GS) (solid blue and solid green curve) and 
for the initial coherent state |acat)lCc3t) (dashed blue and 
green curves and solid red and orange curves). We also 
show the undriven situation (pink and cyan) and the 
static (brown dashed curve). The behavior of (a^a)/j 
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Figure 39: The scaled mean photon number on a resonance 
(jj = LJo — 1.0 for A = 0.3 and A = 1.0 as a function of the 
driving velocity. We compare (a^a)(i/, S,p)/j for two different 
initial states: |GS) and |acs,t)|CcHt) for the driven and undriven 
situation (with tf = 2n /S,/, in both cases). 



for the two different initial state |GS) and |acst)ICc3t) is 
qualitatively similar. Namely that if the equilibrium crit- 
ical paraboloid is not encircled A < Ac = ^wwo /2 then 
the mean photon number is zero. But if we encircle it 
(A > Ac = ^wwo '/2) the mean photon number becomes 
macroscopic. There is no critical driving velocity in the 
driven case. We remind that for the stationary state 
of the rotated Dicke model (see black curve in Fig 



i:(rot) 



39) 



4A^ 



- Wq. 



Next in Fig. 
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we consider {{k^a))T/j as a function 
of 5^ averaged over a time span of = 150 • 27:/ 5^ at 
resonance = wq = 1.0 for the protocol which encircles 
the equilibrium critical paraboloid A — 1.0. We show 
{{a,^k))T/j for the two initial states |GS) and lac^J ICc^t)- 
Also we compare the driven and undriven situations. 

To summarize, we see that the dynamics of the mean 
photon number for the initial ground state |GS) and for 
the initial coherent state lac^J \ Ccst) ^-re very similar. Evo- 
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aj = 1.0, 1^0 = 1.0, A = 1.0, tj = 150-27r/<5^ 
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Figure 40: The time averaged scaled mean photon number 
on resonance cj = wq = 1-0 for A = 1.0 as a function of 
the driving velocity. We plot {a.^a}{tf,5^)/j for two different 
initial states: \GS} and | ac^t ) I Ccst ) for the driven and undriven 
situations (with tf = 150 • 2n/5^ in both cases). 



Figure 41: Time evolution of the parity operator when the 
dynamics is started from the ground state of the unrotated 
Dicke model |GS). The inset plot shows the dependence on 
the atom-file coupling strength A for different pseudo-spin 
length j. 



lution starting from these initial states for a finite number 
of atoms has a tendency to develop two new sub-phases in 
the super radiant phase. Similarly, evolution from both 
states develops a new dynamical critical coupling A^'^^"'' 
or critical driving velocity S'^^"'^ in the thermodynamic 
limit. This phenomenon is traced back to the interplay 
of geometric and dynamical phases in dynamics. Finally, 
we remark that the non-equilibrium phase diagram for 
the initial coherent state |Q:c,t)ICcst) is the same as for 
the initial ground state |GS). 



4- Parity 11 

Here we illustrate the time evolution of the parity oper- 
ator {ip{t)\Il\ip{t)). The result, for a situation when the 
initial state is the ground state of the unrotated Dicke 
model, |GS), while the evolution is given by HRD(i), is 
shown in Fig. |41| The calculation was done numerically 
according to the Chebyshev scheme for a finite number 
of atoms j. 

One can clearly see that the parity is always constant 
in time. More interesting is that there is no transition 
to zero, at a critical coupling. Apparently the system 
always remains in the state with the same parity equal 
to +1. 



V. CONCLUSION 

In conclusion, we studied the dynamics of the rota- 
tionally driven Dicke model beyond the rotating wave 



approximation. We calculated the time evolution of the 
mean photon number and the expectation values of the 
parity operator of the system for different initial states: 
(i) stationary Dicke coherent state |acat)ICcat)i (ii) the 
stationary coherent state |ac2)ICc2) of the rotated Dicke 
model HrdI^), (iii) Fock state \n = 0)\j,m = —j) and 
(iv) the ground state |GS) of the Dicke model. In or- 
der to understand the influence of the geometric phase 
on the non-equilibrium dynamics we compared the time 
evolution for the rotationally driven Dicke model Hrd (t) 
and the evolution of the Dicke model without driving. 

Here is the summary of our findings for different phys- 
ically relevant initial states, (i) For the stationary Dicke 
initial state (stationary state of the usual Dicke model 
Hd) we observe in the TDL a reentrant meta stable phase 
in the super-radiant phase. This is clearly manifested 
in the time averaged scaled mean photon number. We 
observe that the system developed a dynamical critical 
coupling strength A^'^'*'"^ and dynamical critical driving 
velocity (J^*^"''- Dynamical quantum phase transition is 
identified as a sudden change in the dynamical behav- 
ior of observables as a function of the driving parameters 
in real time. We suggest that this critical behavior re- 
sults form a competition of the geometric phase and the 
dynamical phase which can give rise to a resonance phe- 
nomenon described in Ref. [14 . By looking at the time 
averaged scaled mean photon number we constructed a 
non-equilibrium phase diagram. The time evolution of 
the parity operator reveals a similar phenomenon of ap- 
pearance of a meta stable phase and provides a comple- 
mentary information about a new metastable phase. Fur- 
ther, we observe that this type of non-equilibrium driving 
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allows one to probe the equilibrium quantum criticali- 
ties "from a distance" in parameter space, namely by 
encircling a quantum critical surface in parameter space 
without actual crossing it. This may provide a useful 
experimental hint since the quantum state will not be 
destroyed by crossing the quantum critical point. 

(ii) For the stationary circle initial state |ac2)ICc2) 
(stationary coherent state of the rotated Dicke model 
Hrd) we found a shift in the critical coupling by the 
amount given by the applied rotation velocity A^"^"*^ = 
^Ju}(wQ + (5^) /2. This defines a critical rotation driving 

velocity (5^"^°* ■* = ^ — ujq. Which can otherwise be under- 
stood by making a transformation into a rotating frame 
and considering the effective Hamiltonian Hrqt in this 
basis. The expectation value of the parity operator shows 
a constant time evolution and goes to zero at the critical 

1. -v (rot) 

couphng Ac . 

(iii) For the Fock state with no photons and lowest 
spin projection \n = 0)|j, —j) we obtain, for a finite num- 
ber of atoms, a mean photon number which is different 
from zero. We observe the same rotational critical cou- 
pling Ac "^'^ = ■\/aj(a;o + 5^) /2 in this case and there- 
fore (5^™''' ~ ^ Wo, the same as for the initial state 

|ac2)|Cc2)- The mean field solution does not evolve in 
time and produces trivial results if we choose the param- 
eters of the initial conditions as ao = ao = Co = Co = 0. 
However, calculations for a very small initial value of 
the parameters (which can be interpreted as quantum 



fluctuations) ao = ckq = Co = Co = 10"*^ reveal a 
non zero mean photon number in the TDL using the 
mean field equations. For the time averaged quanti- 
ties we also observe the same rotational critical coupling 
aI™'^ = ^/uj(y!vi + 5^)' 12 and (5^™*'' = ^ - wq, as be- 
fore. The parity is always constant and equals one for a 
finite number of two-level atoms which also holds in the 
TDL. On the other hand, for small non-zero initial con- 
ditions we observe an oscillatory behavior of the parity 
expectation value. 

(iv) For the initial state given by the ground state 
|GS) (which is computed numerically by diagonalizing 
the Dicke Hamiltonian Hd) we observe the same dynam- 
ical behavior as for the stationary Dicke state \ac^i)\C,Cst) ■ 
This shows that these two states become equivalent in 
the TDL. 

Summarizing, non-equilibrium dynamics of the Dicke 
model shows several universal properties independent on 
the precise nature of the initial state and exhibits inter- 
esting non-equilibrium quantum critical behavior which 
can be traced back to a competition of geometric and 
dynamical phases. 
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